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We present high resolution single crystal time-of-flight neutron scattering measurements on the 
candidate quantum spin liquid pyrochlore Tb2Ti2 07 at low temperature and in a magnetic field. At 
^70 mK and in zero field, Tb2Ti207 reveals diffuse magnetic elastic scattering at positions 
in reciprocal space, consistent with short-range correlated regions based on a two-in two-out spin 
ice configuration on a doubled conventional unit cell. This elastic scattering is separated from very 
low energy magnetic inelastic scattering by an energy gap of 0.05 meV, which is clearly resolved 
here. The elastic signal disappears under the application of small magnetic fields and upon elevating 
temperature. Pinch-point-like elastic diffuse scattering is observed near (1,1,1) and (0,0,2) in zero 
field at ^70 mK, in agreement with Fennell et al. pQ, supporting the quantum spin ice interpretation 
of Tb 2 Ti 2 7 . 



PACS numbers: 75.25.-j 75.10.Kt 75.40.Gb 75.40.-s 

Geometrically frustrated magnets have been intensely 
studied over the course of the last decade due to the 
rich variety of unconventional ground states they display 
[2]. In particular, the magnetic rare-earth titanate py- 
rochlore oxides with the chemical formula R2Ti2C>7 have 
received considerable attention [3]. In this structure, 
trivalent R 3+ rare-earth metal ions occupy a network of 
corner-sharing tetrahedra known as the pyrochlore lat- 
tice, which is the prototypical example of geometric frus- 
tration in three dimensions. The interplay between ex- 
change and dipolar interactions with crystal-field (CF) 
induced anisotropy on the underlying pyrochlore lattice 
leads to a variety of exotic phenomena. These include 
ground state selection by order- by-disorder [4j[5], as we U 
as classical [6] and quantum spin ice [7] physics. 

Tb2Ti207 has attracted much interest as a potential 
experimental realization of a quantum spin liquid [8], 
based on its lack of long range order (LRO) down to 
at least 50 mK [9j [10| despite a Curie Weiss tempera- 
ture of — 14 K [IT]. The Tb 3+ ions in Tb 2 Ti 2 7 have 
an Ising CF doublet ground state that leads to <111> 
easy-axis anisotropy as in the canonical spin ice materials 
Ho2Ti207 and Dy2Ti2C>7. In contrast to these spin ices 
however, the net exchange in Tb2Ti207 is antiferromag- 
netic (AF). Since the combination of Ising-like <111> 
anisotropy and AF interactions on the pyrochlore lattice 
is unfrustrated [12], Tb2Ti207 is naively expected to dis- 
play a unique LRO ground state. 

Recently, two theoretical pictures have sought to ex- 
plain its enigmatic ground state. One of these proposed 
a "quantum spin ice" (QSI) scenario for the disordered 



ground state in Tb2Ti 2 07 p~3j[T4]. This scenario invokes 
virtual quantum excitations between the CF ground state 
doublet and the first excited doublet (A ~15 K) that 
reposition Tb2Ti207 into the spin ice regime. However, 
Jahn- Teller physics and a non-magnetic singlet ground 
state have also been proposed as the cause underlying its 
failure to order [T5] . 

In this Letter, we report new neutron scattering mea- 
surements on single crystalline Tb2Ti207 that elucidate 
the nature of its ground state. Our measurements were 
performed in zero field and in magnetic fields applied 
along two high-symmetry directions in the pyrochlore lat- 
tice. The zero field measurements reveal short-range AF 
spin ice correlations extending over roughly two conven- 
tional unit cells, characterized by elastic diffuse scatter- 
ing around (|,|,|) positions in reciprocal space. The 
previously reported checkerboard-pattern of diffuse scat- 
tering [T6l[T7] is now resolved as low-energy inelastic scat- 
tering for 0.05 meV < E < 0.6 meV, while the higher- 
energy inelastic scattering regime for 0.8 meV < E < 
2.0 meV is dominated by CF excitations, in agreement 
with earlier studies [Til HE] • A model of the elastic scat- 
tering intensity involving 128 spins in a doubled conven- 
tional unit cell is presented. A fit of this model to the 
observed peak intensities suggests that the Tb 3+ mo- 
ments form an AF ordered spin ice configuration with 
spins tilted ~12° from their local <111> axes. Finally, 
pinch-point like diffuse elastic scattering is observed near 
(0,0,2) in zero field. These results support the QSI pic- 
ture of Tb 2 Ti 2 7 for /i H = T and T-70 mK, and 
suggest that Tb2Ti207 may well form a long range or- 
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FIG. 1. Neutron scattering data within the (H,H,L) plane of Tb 2 Ti 2 7 at T = 70 mK are shown for a) -0.1 meV < E < 0.1 
meV, b) 0.1 meV < E < 0.6 meV and c) 1.0 < E < 1.8 meV. Panel d) shows a plot of intensity vs energy transfer at the 
(§>§>§) position for fields /ioH=0,2,3, and 4T. The integration range is H= [0.2,0.8] meV, L= [1.2, 1.8] meV. The inset expands 
the intensity scale. All data shown were corrected for detector efficiency, and an empty can background was subtracted. The 
error bars are ±lcr. 



dered (LRO) equilibrium state based on this structure 
at lower temperatures, though this might require more 
pristine samples than are currently available. 

The single crystal sample of Tb2Ti207 used for both 
neutron scattering measurements was grown using the 
optical floating zone technique at McMaster University 
P~9j [20]. It is the same single crystal used in an ear- 
lier study by Rule et al. [17]. Time-of- flight neutron 
scattering measurements were performed using the disk- 
chopper spectrometer DCS [21] at the NIST Center for 
Neutron Research and the LET spectrometer [22] at the 
ISIS Spallation Neutron Source. For the DCS measure- 
ments, incident neutrons of energy E^=3.27 meV were 
employed, giving an energy resolution of 0.1 meV. The 
sample was carefully aligned with the [1-10] direction ver- 
tical to within 0.5 degrees, such that the [H,H,L] plane 
was coincident with the horizontal scattering plane. In 
the LET experiment, high resolution measurements us- 
ing E^=2.32 meV gave an energy resolution of 0.02 meV. 
At LET, the sample was aligned with the [11-1] direction 
vertical, placing the plane spanned by [H,H,2H] and [H,- 
H,0] within the horizontal scattering plane. Both DCS 
and LET experiments achieved a base temperature of 
~70 mK and maximum magnetic fields of /ioH = 10 T 
and /ioH = 7 T, respectively. 

Figures l(a)-l(c) show reciprocal lattice maps in the 
[H,H,L] plane in zero field, taken at T~ 70 mK on DCS. 
The data were corrected for detector efficiency, and to 
eliminate scattering from the sample environment. Fig- 
ure 1(a) shows data in the elastic channel, integrating 
over -0.1 meV < E < 0.1 meV. These data reveal strong 
diffuse elastic scattering at (|,|,|)-type positions in re- 
ciprocal space. The isotropic width of the scattering im- 
plies a correlation length of £ ~10 A which corresponds 
to a correlated region of diameter 2£ = 20 A or roughly 
two conventional unit cells. Figure 1(b) shows the low- 
lying inelastic scattering for 0.1 meV < E < 0.6 meV. 
Within our energy resolution, this scattering appears to 



be nearly quasi-elastic and forms the distinct checker- 
board pattern observed previously [16] with highest in- 
tensities near (0,0,2) and (2,2,0) in the [H,H,L] plane. 
Interestingly, this feature has boundaries formed by the 
(|,|,|) diffuse elastic scattering. The higher energy in- 
elastic scattering for 1.0 meV < E < 1.8 meV is shown in 
Fig. 1(c). This intensity results from excitations to the 
lowest excited CF doublet states at A ~1.2 meV. 

Figure 1(d) shows the intensity around the (§,§,§) po- 
sition as a function of energy for selected magnetic fields. 
The zero field data nicely illustrate the three distinct en- 
ergy regimes which are mapped in reciprocal space in 
Figs. 1(a)- 1(c). Under the application of a 2T field, the 
zero field elastic scattering at the (|,|,|) and related po- 
sitions drops to one seventh of the initial intensity, as the 
elastic magnetic scattering has been eliminated, leaving 
only nuclear incoherent scattering at the (§,§,§) posi- 
tion. The low-lying inelastic scattering displays a clear 
peak at ~0.2 meV at 2T. This feature moves to higher 
energies and fades in intensity in larger magnetic fields. 
The broad, higher energy inelastic scattering from the CF 
excitations also splits in the presence of a [1-10] magnetic 
field, forming spin wave bands which were previously re- 
ported [T7] , 

The new elastic diffuse scattering at (|,|,|) and re- 
lated wavevectors can be well described by a model based 
on a "2-in 2-out" spin ice arrangement which extends over 
two conventional unit cells in all three dimensions, lead- 
ing to a supercell composed of 8 conventional unit cells. 
This model takes into account a) the observed phase shift 
of 7r between neighboring unit cells giving rise to scatter- 
ing at (§,§,§) positions and b) the observed width of 
the diffuse scattering. The spin arrangement over two 
neighboring cells in the xy-plane is shown in panel c) of 
Fig. 2. Note that the black spins are reversed in their 
neighboring cell (grey spins) along the x- or y-directions. 

We calculate the neutron scattering intensity S(Q) 
which is proportional to the square of the transverse 
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FIG. 2. Comparison of the measured elastic diffuse scattering 
in the (H,H,L) plane of Tb 2 Ti 2 7 at 70 mK and /i H = 
(panel a)) to the calculated S(Q) discussed in the text (panel 
b)). The intensity scale in the calculation is normalized to the 
largest intensity at (§,§,§)• The spin arrangement between 
neighboring conventional unit cells is shown as a projection 
onto the xy-plane in panel c). The ordered 2-in 2-out spin 
configuration in a single unit cell is shown in panel d), and 
this pattern is reversed in the neighboring cells (black vs grey 
arrows in panel c)) to form "< ~, |, - > ordered spin ice". 
The spins are tilted from their local <111> axes by ~12°. 



component of the magnetic structure factor |M^(Q)| 2 = 
|QxM(Q)xQ| 2 withM(Q) = J2j /j(Q).m j .e^-, where 
j runs over the 128 spins in the supercell. S(Q) is then 
averaged over all three spin ice domains. This calcula- 
tion includes the magnetic form factor fj(Q) for Tb 3+ 
and allows for direct comparison with experimental data 
[23] . In order to determine the deviation of the local 
structure from an ideal 2-in 2-out scenario, we fit the 
calculated intensity to that at the seven (^,^,^)-related 
positions observed in the experiment, letting the cant- 
ing angles from <111> vary. The best agreement with 
the experimental data (Fig. 2 a)) is shown in Fig. 2 b). 
The calculated magnetic structure factor reproduces the 
net intensities of the diffuse (^,^,^)-like peaks very well. 
The elongation of the calculated diffuse scattering is a 
result of the finite size of the assumed spin arrangement 
(128 spins) for which we calculated the Fourier trans- 
form. The best-fit spin configuration has all spins canted 
by ~12° from their local <111> axis with a reduction 
of the magnetic moment along the local z-axis. Our fit- 
ting results are not very sensitive to the canting of the 
other two spin components. We show the resulting spin 
arrangement in Fig. 2 d) for one conventional unit cell. 
Tb2Ti207 at 70mK is therefore composed of short range 
ordered (SRO) domains of an AF < \->\-> \ > ordered 
spin ice. 

The SRO AF spin-ice state of Tb 2 Ti 2 7 at -70 mK 
and zero field shows Tb2Ti207 to have a strong tendency 
to order. An actual LRO state may require lower tem- 
peratures and equilibrium conditions. Sensitivity to weak 



disorder, as is known to occur in single crystals of QSI 
Yb2Ti 2 07 [24], may also be relevant. We note that the 
spin canting angle from the local <111> axis is similar 
to the — 13° spin cant observed in the ferromagnetic "or- 
dered spin ice" phase in Tb2Sn207 [25 . We further note 
a q=(|,|,^) ordering wavevector is known to be selected 
by the combination of isotropic near neighbour exchange 
and dipolar interactions, as occurs in Gd2Ti207 [26] . 
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FIG. 3. High energy resolution data from LET, with the 60K, 
0T data set subtracted, a) Elastic scattering in the [K,K,-K] 
vs. [|H,|H,H] plane at /irjH=0, integrated from E=[-0.02, 
0.02] meV, and [H,-H,0] = [-0.6,-0.4] r.l.u.. In the scattering 
geometry used in the LET experiment, the [K,K,-K] direction 
is parallel to the vertical (field) direction. The diffuse scatter- 
ing at 66mK, 0T shown here corresponds to (-§,§,-§)? which 
arises at [|H,|H,H] = -0.33 r.l.u., [H,-H,0] = -0.5 r.l.u., and 
[K,K,-K] = 0.166 r.l.u.. This diffuse elastic scattering disap- 
pears upon warming to 1.5K. b) Field and temperature de- 
pendence of the diffuse elastic scattering at (- - , \ - \ ) as shown 
in panel a), integrating over [|H,|H,H] = [0,0.3] r.l.u.. c) In- 
elastic scattering integrated over [|H,|H,H] = [-1.5,1.0] r.l.u. 
and [K,K,-K] = [-0.5, 0.5] r.l.u.. A 0.05 meV gap opens up 
below 1.5K, which was not resolved previously, d) Field and 
temperature dependence of the inelastic scattering as shown 
in panel c), integrated over the full [H,-H,0] range. 

Figure 3 shows the field- and temperature-dependence 
of the scattering near the (~\,\-\) position as measured 
on LET with an incident neutron energy of 2.32 meV. 
A high-temperature, paramagnetic data set at T = 60 
K has been subtracted from all data sets in this figure. 
Fig. 3 a) shows elastic scattering, binned between -0.02 
meV < E < 0.02 meV. The diffuse elastic scattering is 
of the same (|,^,|)-type that was measured on DCS, 
although it occurs here within a different scattering plane 
than was shown in Figs. 1 and 2. Fig. 3 b) shows 
cuts through this (-§,§,-§) position for different fields 
and two temperatures. The (-|,|,-|) elastic scattering 
can be completely removed by an applied field of 0.25 T 
along [11-1], confirming that it is magnetic and of the 
same origin as the elastic scattering shown in Figs. 1 
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and 2. The elastic scattering also vanishes by T = 1.5 K, 
indicating that it is likely related to the broad hump in 
specific heat observed at 0.4 K in many samples [TTJ|27]- 

GE|. 

Figure 3 c) shows energy vs. wavevector slices of these 
same data at T = 1.5 K and 0.07 K with /iqH = 0, as 
well as T = 0.07 K and /i H = 0.25 T and 7 T. These 
new, high energy resolution data clearly show the "quasi- 
elastic" spectrum to be suppressed at low energies and 
T = 0.07 K and /i H = T relative to that at T=60 K; 
that is to display a gap of ~ 0.05 meV (i.e. 0.57 K). The 
gap is itself very sensitive to the magnetic field, expand- 
ing to - 0.08 meV for /i H = 0.25 T. Fig. 3 d) shows 
the energy dependence of this scattering (integrated over 
the full range measured in the two perpendicular Q di- 
rections), relative to that at T = 60 K. The clear pres- 
ence of the ~ 0.05 meV gap separating the low energy 
inelastic scattering (0.05 meV to 0.6 meV) from the elas- 
tic scattering due to the SRO AF spin ice scattering at 
(|,|,|) wave vectors, confirms that three distinct regimes 
of magnetic scattering exist at low temperatures: elastic, 
low energy inelastic, and "high" energy inelastic scatter- 
ing due to the excited CF doublet at energies of ~ 1 meV 
(see Fig. 1 d). Our observation of elastic magnetic scat- 
tering below 0.05 meV in zero field is inconsistent with 
the non-magnetic, singlet ground state scenario previ- 
ously suggested [15] . 
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FIG. 4. Pinch-point like scattering in the ground state of 
Tb 2 Ti 2 7 observed on DCS at - 70 mK for /x H = 0T in 
panel a). Panel b) shows diffuse scattering along the <0,0,L> 
direction under application of a /ioH = 2T field, which might 
originate from anisotropic exchange. Both data sets show 
elastic scattering for -0.1 meV < E < 0.1 meV and are dis- 
played on a logarithmic intensity scale (compared to a linear 
scale in Fig. 1). 

Pinch-point diffuse scattering characteristic of 
Coulombic spin ice correlations [31- 33 has been of great 
interest in the classical spin ice materials Dy2Ti20y and 
Ho2Ti207 and recently in Tb2Ti2C>7 pQ. Figure 4 shows 



these correlations, most notably around (0,0,2), as they 
appear in the elastic scattering at T = 0.07 K and (a) 
zero field, and with a (b) 2 T magnetic field applied 
along [1-10]. Our pinch point diffuse scattering near 
(0,0,2) is qualitatively similar to that observed in Fig. 
2a) of Fennell et al. pQ. As shown in Fig. 4, we observe 
"rods" of scattering along the <0,0,L> and <H,H,H> 
directions, which appear to be modulated by pinch 
points near the zone centers at (1,1,1), (1,1,3) and near 
(0,0,2). In a 2 T field applied along [1-10], all elastic 
scattering at (|,|,|) positions has disappeared and only 
broad diffuse scattering along the <0,0,L> direction 
remains. The origin of this rod-like elastic diffuse 
scattering is not understood, but one can speculate that 
it originates from anisotropic exchange in Tb2Ti207, as 
is known to describe the microscopic spin Hamiltonian 
in Yb 2 Ti 2 7 [7 and Er 2 Ti 2 7 @]. 

The ground state of Tb2Ti2 07 has by now been studied 
for more than a decade by neutron scattering techniques. 
It is surprising that the (^,|,|) elastic SRO scattering, 
presaging a transition to LRO, has not previously been 
observed. This is likely due to the low onset temperature 
of this scattering: for example Rule et al. [17] previously 
studied the same single crystal with the same instru- 
mental conditions at DCS, but at temperatures of 400 
mK and above, without seeing this scattering. In fact, 
triple- axis measurements by Yasui et al. [34 do show 
evidence for a peak at (|,|,§) at T = 50 mK (Figure 
2 in [34 j), however this feature was not pursued further 
in subsequent studies. Diffraction experiments, such as 
that by Gardner et al. [16], would be problematic due 
to integration over elastic and low energy inelastic scat- 
tering, which we show in Fig. 1 to have very different 
Q-dependencies. 

In conclusion, new neutron scattering measurements 
on Tb2Ti207 at T = 0.07 K in zero field have revealed 
elastic diffuse scattering at (|,|,|) positions characteris- 
tic of short-range AF spin ice correlations extending over 
roughly two conventional pyrochlore unit cells. This elas- 
tic scattering is separated by a gap of ~ 0.05 meV from 
low lying inelastic scattering, and can be quantitatively 
described based on an ordered 2-in 2-out local spin ice 
structure with a spin canting angle of ~ 12°. This devel- 
opment of AF spin ice correlations in Tb2Ti207 is char- 
acteristic of a strong tendency to form the corresponding 
LRO state which does not occur either because the tem- 
perature is not sufficiently low, or due to weak disorder in 
the samples, as has been recently characterized in other 
exotic pyrochlore magnets [24], or both. 
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